Theoretical study of a cold atom beam splitter 



O 
O 

(N 

wo: 
< 

(N 
> 

m 



O ■ 
\o ■ 
o ■ 

^ : 

9 L,: 

^— > ' 
G ■ 
cd 

3 : 
cr. 



Naceur Gaaloul, 1,2 Annick Suzor-Weiner, 1 Laurence Pruvost, 3 Mourad Telmini, 2 and Eric Charron 1 

1 Laboratoire de Photophysique Moleculaire, CNRS, Batiment 210, 
Universite Paris-Sud 11, 91405 Orsay cedex, France. 
a Laboratoire de Spectroscopie Atomique, Moleculaire et Applications, 
Department of Physics, Faculty of Sciences of Tunis, 
University Tunis El Manar, 2092 Tunis, Tunisia. 
3 Laboratoire Aime Cotton, CNRS, Batiment 505, 
Universite Paris-Sud 11, 91405 Orsay cedex, France. 
(Dated: February 1, 2008) 

A theoretical model is presented for the study of the dynamics of a cold atomic cloud falling 
in the gravity field in the presence of two crossing dipole guides. The cloud is splitted between 
the two branches of this laser guide, and we compare experimental measurements of the splitting 
efficiency with semi-classical simulations. We then explore the possibilities of optimization of this 
beam splitter. Our numerical study also gives access to detailed information, such as the atom 
temperature after the splitting. 

PACS numbers: 03.75.Be ; 32.80.Pj ; 32.80.-t ; 39.25.+k 



I. INTRODUCTION 

The manipulation of cold atoms with optical fields is 
a very promising technique which is rapidly developping 
in the context of atom optics Q. Its applications range 
from laser cooling and trapping 0, 0, S @ to coherent 
atom tr ansp ort [3 0- IToUllj and matter wave interfer- 
ometry [12], [l3| ■ Optical fields have also been proposed 
as an interesting tool to control the dynamics of internal 
and motional states of cold atoms for quantum informa- 
tion processing 0, 0, 0, 0, 0, ^1 , and recent exper- 
imental studies have demonstrated a very pr omising first 
implementation with optical lattices |20l |21| . 

For these types of applications, an effective way of 
guiding or transporting the atoms while keeping their 
coherence is required. For atom interferometry it is also 
necessary to separate the atomic wavefunction between 
the arms of an interferometer. Several experimental con- 
figurations have thus been explored for the implementa- 
tion of an atom beam splitter with optical 0, E3 and 
magnetic 013 HHH field potentials. 

Recently, various configurations of atom beam split- 
ters have been implemented and studied with coherent 
sources of atoms [2(| [53, [28( . Indeed Bose- Einstein con- 
densates are the best candidates for the implementation 
of cold atom interferometers because of their intrinsic co- 
herence. However, for certain applications such as atomic 
clocks, gyrometers and gravimeters, thermal ensembles 
of cold atoms are still preferred because they allow for 
the preparation of a larger number of atoms at a higher 
repetition rate Furthermore, recent experiments [|J 
show that high-contrast interference patterns with fringe 
visibility greater than 90% can be obtained with thermal 
sources of atoms. 

In the case of dilute condensates, the dynamics is eas- 
ily followed by solving the one-, two- or three-dimensional 
time-dependent Gross-Pitaevskii equation [2^, |3(j ■ With 
a thermal source of atoms, the theoretical description is 



more complex in the sense that many transverse modes 
of the guide are populated initially in an incoherent way. 
Implementing a cold atom interferometer with this type 
of atom source is therefore not trivial and requires a de- 
tailed theoretical investigation. In addition, new oppor- 
tunities are now investigated with spatial light modula- 
tors [3lJ for the implementation of atom optics devices 
(guides, splitters, . . . ) for thermal and coherent sources 
of atoms. The development of new theoretical models 
aimed at the description of the dynamics of a guided 
thermal ensemble of cold atoms is therefore needed for 
these applications. 

A few theoretical studies of cold atom beam splitters 
have been published recently in various trapping situa- 
tions 13 . In these approaches it was assumed 
that the atomic wave packet was tightly confined in one 
dimension, and the effect of gravity was neglected. In 
the present article we study the cold atom beam splitter 
implemented in reference |jj , by solving numerically the 
time dependent Schrodinger equation for the atomic mo- 
tion in the presence of the gravity field and with realistic 
trapping potentials. Our aim is to propose a theoretical 
model which can reproduce the main features of this ex- 
periment, and to test its pertinence for more elaborate 
situations. 

A large ensemble of 87 Rb atoms is initially trapped and 
cooled to a temperature To in a magneto-optical trap 
(MOT) localized at the height z = (see Figure Q for 
a schematic view). At time t = 0, the trapping poten- 
tial is switched off, while a vertical far off-resonant laser 
beam, crossing the cloud close to its center, is switched 
on. The dipole interaction creates a potential well of 
depth Uq which traps a significant portion of the atoms 
in the transverse directions x and y. 

The guided atoms then fall due to gravity, with a con- 
fined dynamics in the x and y directions. At time to, a 
second oblique guide is suddenly switched on. The two 
guides cross at the height z c and form an angle 7. A po- 
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FIG. 1: Schematic view of the guiding setup. The initial 
87 Rb cloud is shown at z = 0. When the magneto-optic trap 
is switched off, a vertical laser beam is switched on and the 
cloud, partially trapped by the associated dipole force, falls 
in the gravity field. At time to a second oblique guide is 
switched on. The two guides cross at the height z c and form 
an angle 7. The trapping potential seen by the atoms at two 
different heights is shown in the insets. Finally, the atoms are 
probed 1 cm below the initial position. 

tential well, of depth U\, is induced by the optical dipole 
interaction with the oblique guide. This creates an addi- 
tional path for the motion of the atoms. Depending on 
the various parameters (light intensity, angle 7, temper- 
ature To, ... ) a splitting of the cloud can be observed Q 
between the vertical and oblique branches. 

II. THEORETICAL MODEL 

In order to simplify the numerical treatment of this 
phenomenon, we restrict the dimensions of this study 
to the plane defined by the two guiding beams, i.e. to 
the x and z dimensions only. In addition, we adopt a 
semi-classical approach, where the effect of the gravity is 
treated classically. This approximation is justified by the 
value of the de Broglie wavelength associated with the 
speed of the particles in the z direction, Adb *~ 1 A. 

A. Classical approach 

A two-dimensional classical trajectory {x(t), z(t)} is 
first evaluated by solving Newton's equations of motion 



for an atom initially at the position {xq, zq] with the mo- 
mentum {io,io}- An efficient variable time-step Runge- 
Kutta integrator |3!| is used to solve these equations in 
the total potential Vt(x, 2, t) — V g (x, z, t) + mgz, where 
m denotes the atomic mass and g the gravitational con- 
stant. The guiding potential V g (x, z, t) is given by the 
following sum 

V g (x,z,t) = V (x) + Vi(x,z,t), (1) 

with 

f V (x) = -U e~ 2 * 2 /< 

\Vi(x,z,t) = -Ui u(t- t ) e- 2x ' 2 / w * [ 1 

In the previous expressions wq and w\ denote the waists 
of the vertical and oblique laser beams respectively, and 
u{t — to) stands for the Heaviside step function. The 
following rotated coordinates 

J x' — x cos 7 + (z — z c ) sin 7 ,_, 
[ / = (z — z c ) cos 7 — x sin 7 * ' 

have also been introduced (see Figure^)- The Gaussian 
form of this potential arises from the Gaussian intensity 
profile of the laser beam 

Figure^ shows three typical trajectories with poten- 
tial parameters close to the one chosen in the experiment 
performed in Orsay pj ■ From this graph it is already clear 
that, classically, only very specific initial conditions (such 
as xq = —0.2 mm and zo = xq = Zq — 0) drive the atom 
in the oblique branch. Comparing the kinetic energy of 
the atom in the x' direction at time t — to, E x i, with 
the binding energy of the oblique guide gives the actual 
criterion which decides in favor or against the deviation 
of the atom from its natural vertical fall. This can be 
inferred from Figure [2k : the trajectory induced by the 
initial condition xo = +0.2 mm crosses the oblique guide 
almost perpendicularly at the point labelled 2 on the 
graph, with a maximum kinetic energy E x i . This trajec- 
tory therefore remains almost unaffected by the presence 
of the oblique laser beam. On the contrary, the trajec- 
tory associated with the initial condition xo = —0.2 mm 
meets the oblique guide almost tangentially at the point 
labelled 3, and the atom is deflected from its initial ver- 
tical motion. The intermediate case xo = follows a 
dynamics similar to the initial condition xq = +0.2 mm 
with a very slight deviation from the initial vertical mo- 
tion. It will be shown in section II I II that this simple 
interpretation, in terms of individual trajectories, is not 
valid anymore when the atomic dynamics is treated at 
the quantum level. 

Finally, Figure |2b shows the deviations of these 
trajectories from a simple free fall dynamics in 
the z-direction defined by the classical expres- 
sion Zff(t) = zo — Zot — \gt 2 ■ On the time interval 
^ t ^ 45 ms, required to reach the detection probe 
located at z p = —10 mm, this deviation remains smaller 
than half a percent, and this justifies the approach 
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FIG. 2: (Color online) (a) Classical trajectories of a cold atom 
falling in the gravity field in the presence of the two trapping 
potentials with Uo — 30 ^K, Ui = 10 /j,K, wo = 0.2 mm and 
w\ = 0.3 mm. The three trajectories correspond to an ini- 
tial height zo — and initials positions xo = —0.2 mm (black 
solid line), xq — (red dashed line) and xo = +0.2 mm (blue 
dash-dotted line) with zero initial momentum (io = and 
io = 0). The two guides cross at the height z c — —4 mm with 
an angle 7 = 0.12 rad, and the oblique guide is switched on 
at to = 28.6 ms. At this time, a free fall dynamics would give 
the position Zff(to) = Zq — Zot — \gt'o = —4mm correspond- 
ing exactly to z c . Unless specified, these laser parameters 
remain fixed throughout the paper. The thin vertical and 
oblique dotted lines reveal the geometry of the laser beams, 
(b) Deviation of the trajectories from a free fall dynamics in 
the 2-direction for the same initial conditions as in (a). 



thermal average once &(x, z, t) is known (this averaging 
procedure is explained in section 11*1*11 Eq. I|19l21fl ). The 

two-dimensional Hamiltonian H.2F>(x, z, t) can be written 
as the following sum 



H 2 b (x, z,t) — T x + T z + V g (x, z, t) + mgz , (5) 



where 



q 2m dq 2 



(6) 



denotes the kinetic energy operator along the q- 
coordinate. 

The semi-classical approximation discussed above is 
introduced by replacing the z coordinate in the two- 
dimensional Hamiltonian H.2d(x, z,t) with the simple 
classical parameter Zff(t) = zq — z^t — \gt 2 . The quan- 
tum dynamics of the atomic cloud along x is then de- 
scribed by the time-dependent Schrodinger equation for 
the wavepacket (p(x,t) 



d 

ih—<p(x,t) = Hm(x,t) <p(x,t) , 
ot 



(7) 



with a one-dimensional time-dependent Hamiltonian 
H\T>{x,t) given by 



Hm{x,t) = T X + Vm(x,t) 



(8) 



where Vir>(x,t) = V g (x, Zff(t),t). By adopting this 
one-dimensional approach, the numerical simulation is 
simplified at the cost of replacing the two-dimensional 
potential V g (x,z,t) (equations ("JJ and J""J), which is 
time-independent for t > f , with the time- varying one- 
dimensional potential Vix>(x, t). This potential changes 
slowly during the fall of the atom, and two of its snap- 
shots are shown in the insets of Figure []] 



adopted in the following, where the dynamics along the 
z-dimension is simply treated as a classical free fall. 
We now turn to the description of this semi-classical 
approach. 

B. Semi-classical treatment 

In the limit of dilute gases, the dynamics of the 
cloud can be simulated by solving the time-dependent 
Schrodinger equation along the x and z dimensions for 
the wavepacket ^(x, z, t) describing the external dynam- 
ics of a trapped atom 

d 

ih—^{x, z, t) = H 2 t>(x, z, t) z, t) . (4) 
at 

Since the initial state of the atomic wave packet can be 
described in general by a thermal mixture \32\ , the calcu- 
lation of an observable at time t can be done by a simple 



C. Time-dependent propagation 

We assume the atom to be initially in a well defined 
vibrational level vq of the vertical guide potential Vq (x) 

ip(x,t = 0) =Xv (x), (9) 

and we propagate the translational wavepacket until the 
time tf corresponding to the height of the detection 
probe, using the splitting operator method developed by 
Feit et al [36| 

tp{x, t + St) = c-^mst/h ^ t) _ (10) 

The total Hamiltonian 7iiT>(x, t) is splitted in two parts 
corresponding to the kinetic and potential propagators 

~m 1B st/h _ -if x 6t/2h -iV 1B (x,t)8t/h -if x st/2h 

c — C Ac Ac 

(11) 

to decrease the error to the order (6t) 3 . The kinetic 
propagation is performed in the momentum space, and 
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the potential propagation in the coordinate space. Fast 
Fourier Transformation (FFT) allows rapid passage back 
and forth from one representation to the other at each 
time step. Typical grids extend from x m i n = —1.0 mm 
to x max — 2.0 mm with N — 2 20 grid points, and a time 
step of the order of St ~ 40 /is is used. 

At the end of the propagation the wavefunction 
ip(x,tf) is analyzed to determine the efficiency of the 
beam splitter and to extract detailed information on the 
state of the atom in each branch of the laser guide. 



D. Initial trapping of the atomic cloud 

We assume the initial atomic cloud to be in a ther- 
mal state at temperature To described by the usual 
Maxwcll-Boltzmann phase-space probability distribution 
W(x, z,x, z) defined by the following four-dimensional 
product 

W{x,z,x,z) = W Q (x)xW Q (z)xWp(x)xWp(z) , (12) 

where Wq{q) and Wp(q) (q = x,z) are the position and 
momentum distributions 



W Q (q) 
W P {q) 



1 



2tt k B T 



-q 2 /2o 



-mg 2 /2fe B T 



(13) 



In this expression <7o characterizes the size of the cloud 
and k B is the Boltzmann constant. When the vertical 
guide is suddenly switched on at time t = 0, only a frac- 
tion of the cold atoms are trapped in the dipole potential 
created by the laser beam intensity profile. 

The total trapping probability Ptrap can be calculated 
from the position and momentum distributions, assuming 
that for a given position x the atom is trapped if its 
kinetic energy E x along this direction is lower than the 
binding energy —Vq(x) of the potential. The integral over 
x leads to the following rapidly convergent expansion 



n=0 



Ptra P X) ( fco7 



k B T J (2n + 1) n\ ' 



where 
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[wl + {An + 2)al] 



1/2 



Uo 



irk B T 



(14) 



(15) 



This expression is the one-dimensional analogue of the 
well-known two dimensional probability given for exam- 
ple in [U El El 

It should be noticed here that the 
trapping probability only depends on the two following 
dimensionless ratios oq/wq and Uo/k B To- This is a signa- 
ture of the fact that the trapping probability can be ex- 
pressed in a phase-space diagram as the overlap between 
the atomic cloud distribution and the trapping condi- 
tion mm. 
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FIG. 3: (Color online) One-dimensional (black solid line) 
and two-dimensional (red dashed line) trapping probability 
as a function of the dimensionless ratio £/o/&bTo (logarithmic 
scale) for ao/wo = 1.5. The blue point located at the abscissa 
Uo/k-BTo — 1.3 with the error bar is an experimental mea- 
surement extracted from reference [3S|. The one-dimensional 
trapping probabilities for ao/wo = 0.5 and 2.5 are also shown 
as blue and green dotted lines with the labels (1) and (2) 
respectively. All other parameters are as in Figure 



The variation of the trapping probability with the ra- 
tio Uo/ksTo is depicted in Figure El For a given C/q, 
the trapping probability changes very slowly with the 
temperature of the atomic cloud. In a real experiment, 
trapping occurs along both x and y, and P t j ap is there- 
fore shown in Figure[3]to compare with the measurement 
performed in [3g for the ratio C/o/A:bTo = 1.3. Finally, as 
one might intuitively guess, the trapping probability in- 
creases significantly when the size of the atomic cloud <7o 
decreases compared to the laser waist wq characterizing 
the size of the trapping potential. 

Using arguments based on energy conservation, the 
probability for an atom to be trapped in a well defined 
initial vibrational state vq of total energy eo can also be 
calculated using 



+lo 



W Q (x)W E (£ -V (x))dx, (16) 



once the density of states p(eo) in the potential Vq{x) 
is known. The positions x = ±Iq, corresponding 
to 2 Zq = Wq ln(— Uq/eo), are the left and right turning 
points of the level vq, and 



W E (E) = 



e - E I k B T 



\J-kE k B T Q 



(17) 



is simply the analogue of Wp(q), written in terms of ki- 
netic energy E = Sq — Vq(x). The state-dependent prob- 
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FIG. 4: (Color online) Initial population of the various vibra- 
tional levels P(vo) normalized with respect to the population 
of the ground state P(0) as a function of the ratio of their 
energy £o to Uo. The black solid, red dashed and blue dash- 
dotted lines correspond to go /wo = 1.5, 0.5 and 2.5 respec- 
tively. All other parameters are as in Figure [5] The energies 
of the levels v = 0, 2000, 4000, 6000, 8000, 10000 and 11704 
are indicated by the thin vertical dotted lines. 
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FIG. 5: (Color online) Contour plot depicting the time evo- 
lution of the envelope of |<p(:r,£)| 2 on a logarithmic scale as 
a function of x and of z(t) — zo — zot — ^ffi 2 - The atom is 
falling in the gravity field in the presence of the two trapping 
potentials with Uo = 30 /jK, U\ = 10 fiK, wo = 0.2 mm and 
wi — 0.3 mm. The two guides cross at the height z c = —4 mm 
with an angle 7 = 0.12 rad, and the oblique guide is switched 
on at to = 28.6 ms. The vertical and oblique dotted white 
lines reveal the directions of propagation of the laser beams. 
The two white arrows point to the location of the maximum 
probability density at the height 2 = —10 mm. The initial 
level is here vq = 6000. 



ability P{vq) only depends on the three following dimen- 
sionless parameters: Co/wo, Uo/ksTo and Sq/Uq. It fi- 
nally satisfies the relation 

Ptrap = X) P M- (18) 

Figure ^ shows the initial distribution of vibrational 
levels for various ratios ao/wQ. The lowest levels domi- 
nate the distribution, and this is particularly true when 
the size of the atomic cloud cto is smaller than the size of 
the trapping potential wq. Just as with the total trapping 
probability Ptrap shown in Figure [21 P{vq) changes very 
slowly with the temperature. This variation is therefore 
not shown here. 



III. NUMERICAL RESULTS 

A typical quantum dynamics can be seen Figure |S] 
which shows the time evolution of the initial level 
i>0 = 6000 as a function of x and z for the initial con- 
dition zq — zq = 0. This initial state is a stationnary 
state of the vertical guide, and it does not evolve in time 
until it reaches the height z c = — 4 mm where the two 
dipole guides cross. Afterwards, a wavepacket is created, 
and this one evolves inside two main branches, indicated 
by the white arrows. The oblique "trajectory" is guided 
by the oblique laser beam represented by the thin white 



oblique dotted line, while in the vertical branch an oscil- 
lating wavepacket is evolving. The picture obtained here 
is not as simple as the classical trajectories shown in Fig- 
ure [21 since this quantum state cannot be represented by 
a single trajectory. The consequence is that, depend- 
ing on the initial conditions, the atomic wavepacket can 
be delocalized in the two guides simultaneously. A single 
initial quantum state can therefore split coherently along 
two paths separated by macroscopic distances. This ef- 
fect, which is quantum by nature, might open some inter- 
esting perspectives for atom interferometry experiments 
with laser guides. 

The probability Pr(vq, Zq, Zo) of finding the atom in 
the right wing potential well and thus in the oblique guide 
at the height of the detection probe (z p — —10 mm) is 
evaluated for each initial trap state vq by 

Pr(vq,zo,z ) = \(p Zo , i() (x,tf)\ 2 dx. (19) 

/ Oblique 
J Guide 

The wavepacket Lp ZQt z {x,t) is labelled here by the in- 
dexes (zq, Zo) indicating the initial conditions of the sim- 
ulation. The probability Pr(vo,0,0) is shown as a black 
solid line in the upper part of Figure as a function of Vq 
for an atom initially at rest (io = 0) at the height zq = 0. 
Clearly, the lowest energy states (v < 4000) are not devi- 
ated by the oblique guide, and they simply fall vertically. 
The explanation for this effect is simple: their energy is 
too small for them to be trapped in the present oblique 
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guide of depth U\ = 10 fj,K (see for instance the lower in- 
set of Figure representing the guiding potential around 
z~ —8 mm). The eigenstates of vibrational quantum 
number higher than vo ~ 5800 are the only states of total 
energy Sq — 10/uK. In an energy-based first approxima- 
tion, all states with vo ^ 5800 should therefore remain 
unaffected by the beam splitter. In reality, all vibrational 
states experience a quickly varying potential in the vicin- 
ity of z c . They are therefore subjected to non-adiabatic 
transitions to higher or lower excited states which may 
or may not be captured in the oblique guide. This non- 
adiabatic effect is especially important for the highest 
falling speeds, and thus for the initial conditions zq > 
and zo ^ 0. The upper part of Figure therefore shows 
the same probability but incoherently averaged over the 
initial distributions of zq and io 
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In this figure, we see that the states of vibrational quan- 
tum numbers 3000 ^ vq ^ 5800 already have a significant 
probability of splitting. A comparison of the averaged 
probability distribution (Pr) (vq) (red dashed curve in the 
upper part of Figure^! with the probability Pr(vq, zq, io) 
obtained for a single initial condition (zq — zq = 0, black 
solid curve in the upper part of Figure EJ) shows that 
even if the averaging procedure modifies significantly the 
probability distribution, a qualitatively correct descrip- 
tion is already obtained by a single calculation with the 
atom initially at rest, with zq = zq = 0. 

The initial position and momentum distributions of 
the atoms along x is taken into account by balancing the 
averaged probability (Pr)(vq) with the initial trapping 
probability P(vq). The splitting efficiency of the state vq 
is therefore written as 



P s (y ) =PM x (P R )(v ) 



(21) 



This state-dependent splitting probability P s (vq) is 
shown in the lower part of Figure as a function of vq for 
various sizes of the atomic cloud <7o at fixed temperature 
To = 14 /iK. An interesting tendency can be noticed in 
this figure : larger atomic clouds, since they favor the 
initial trapping of higher vibrational levels in the verti- 
cal guide (see Figure 0] for instance), have a higher total 
splitting efficiency, and present a distribution of levels 
clearly shifted to higher energies. As a consequence, the 
average value (vo) of the trapped states is (vq) ~ 7750 for 
(To = 0.15 mm and (vq) — 8140 for <7o = 0.45 mm. 

The "total" splitting efficiency P s at temperature To is 
finally evaluated by averaging over the vibrational quan- 
tum numbers vq according to 



(22) 



where Ptrap is the total trapping probability [Eg J14(l ]. 
A unit splitting probability (P s = 1) would indicate that 
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FIG. 6: (Color online) State-dependent splitting efficiency as 
a function of the initial level «o- Upper graph (a) : the prob- 
ability Pr(vo, 0, 0) of finding the atom in the oblique guide at 
the end of the propagation is shown as a black solid line for the 
initial condition zo = io = 0. The same probability averaged 
over zo and io, (Pr)(vo), is shown as a red dashed line. Lower 
graph (b) : The averaged splitting probability P s (vo) scaled 
by the same factor P(0) as in Figure|l|is shown as a function of 
vo for various sizes of the atomic cloud cto and for To = 14 /iK. 
The cloud sizes cro = 0.15, 0.30 and 0.45 mm correspond to 
the solid black, dashed green and dash-dotted blue curves re- 
spectively. These probabilities have been averaged over the 
initial classical conditions chosen for zo and io- The guide 
parameters are as in Figure The two vertical lines indicate 
the average value of «o in the oblique guide for cro = 0.15 mm 
(black solid line), (vq) = 7750, and oo = 0.45mm (blue dash- 
dotted line), (v ) = 8140. 



all trapped atoms are captured by the oblique guide. A 
perfect beam splitter, whose reflection and transmission 
coefficients equal 0.5, corresponds to P s = 0.5. 

The total splitting efficiency P s of the present beam 
splitter setup has been measured recently in Orsay for 
various heights of the crossing point z c [38| . The conclu- 
sion of this experimental study is that, with the parame- 
ters chosen in Figured a maximum splitting efficiency of 
about 10% is observed around z c ~ —6 mm. Some mea- 
surements made with a smaller waist w\ and a higher po- 
tential depth U\ also show that the variation of P s with 
z c is not symmetric with respect to its maximum value 
z c ~ —6 mm. Our numerical study, which gives a max- 
imum splitting probability of about 15% for the height 
z c ~ —5.2 mm (see Figure 01 is therefore in qualitative 
agreement with this experimental measurement. 

The variation of the splitting probability with z c can 
be explained as follows. If the average position of the 
atomic cloud is much lower than the crossing height 
when the oblique guide is switched on (z//(^o) z c ), 
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z c (mm) 

FIG. 7: Total splitting efficiency P s [En.l|22p] as a function of 
the crossing height z c of the laser beams. This probability has 
been averaged over the initial conditions chosen for zo and Zo, 
and summed over vo- The size of the initial atomic cloud is 
<7o == 0.30 mm. All other parameters are as in Figure |S| 



the beam splitter is inefficient, as can be seen in Figured 
for z c 3> —4 mm. On the opposite side of this graph, for 
z c <C —4 mm, the oblique guide is on when the atoms 
reach z c , but they reach this height with a kinetic en- 
ergy which becomes comparable to - or higher than - 
the binding energy of the oblique guide U±. This ex- 
plains why the efficiency of the beam splitter falls to 
when z c <C —4 mm. 

We have also calculated the variation of the total split- 
ting efficiency P s [Ea. (|22[l ] of this beam splitter with one 
of the most crucial parameter : the potential depth of the 
oblique guide U±. This numerical simulation has been 
performed for various ratios of oblique to vertical beam 
waists ivi/wq. The result is shown in FigureJSJwith fixed 
initial conditions zq = zq = 0. One can notice here that 
the splitting efficiency varies monotonically from to its 
maximum value when U\ varies from to 3 Uq. Depend- 
ing on the value of the waist of the oblique guide, a total 
deflection of the beam can be realized (see for instance 
the case w\ = 1.5 wo and U\ = 3 Uq). A completely sym- 
metric splitting is also predicted when U\ ~ 1.1 L7 and 
wi ^ wq. This last prediction is in agreement with the 
experiment |38l |. 

The dash-dotted blue curve of Figure [S] corresponds to 
a ratio of laser waists (wi/wq = 1.5) very close to the 
experimental one [7j. Our semi-classical model repro- 
duces in this case the experimental splitting efficiency of 
44% (see Figure 4-f of reference 0) when Ux/Uq = 0.95. 
A calculation performed for the same ratio of potential 
depths but with z c = —2 mm gives a splitting efficiency 
of 28.3%, again very close to the experimental value of 
29.2% (see Figure 3-b of reference 0). 
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FIG. 8: (Color online) Total splitting efficiency P B [Eg. 11231 1 
for the initial conditions zq — zq = as a function of the 
ratio of the oblique to vertical potential depths Ui/Uo, and 
therefore as a function of the ratio of the beam intensities. For 
this calculation, Uo is fixed (30 fiK) and Ui is varied. The 
crossing height between the two guides is z c = —4 mm. The 
waist of the vertical beam is too = 0.2 mm, and the splitting 
efficiencies calculated for an oblique waist of Wi — 0.1, 0.2 
and 0.3 mm are shown as black solid, red dashed and blue 
dash-dotted lines respectively. All other parameters are as in 
Figure |7| 



Finally, when the oblique guide is deep enough to in- 
duce a significant splitting of the atomic cloud, a higher 
splitting efficiency can always be obtained by increasing 
w\. The results shown in this figure therefore indicate 
that a high degree of control exists in this type of experi- 
mental configuration since the splitting efficiency can be 
modified at will. 

Figure shows the average transverse energy (direc- 
tions x and x') of the atoms in the vertical and in the 
oblique guide after the splitting : (Eq) and (E\). An 
evaluation of the final energy Eq(vq) in the vertical guide 
is first performed for each initial state vq using the ex- 
pression 



E (vo) = 



Vertical 
Guide 



V^o (%, tf) 'Hid (x, t) ip Vo (x, t f ) dx (23) 



This energy is then averaged over all vibrational levels 
E„ P(v )E (v ) 



(Eo) 



(24) 



where 



Vertical 
Guide 



flofatf) <p V0 (x,t f )dx (25) 
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FIG. 9: (Color online) Average energies (Eo) (black solid line) 
and (Ei) (red dashed line) in the vertical and oblique guides 
as a function of Ui/Uq for the initial conditions zq = zq — 0. 
For this calculation, Uo is fixed (30 fiK) and U\ is varied. The 
crossing height between the two guides is z c — —4mm. The 
waist of the vertical and oblique beams are wo = 0.2 mm and 
wi — 0.3 mm. All other parameters are as in Figure [7| The 
total average energy (E) is also shown as a green dash-dotted 
line. 



is the probability of experiencing a simple vertical fall 
when starting in the initial level Vq. In the oblique guide, 
a similar approach is used to calculate the average energy 
(Ei), but the transverse direction is now x' . A rotation of 
the reference frame is therefore in order. For this calcu- 
lation, the wavefunction tp Vo (x,tf) and the Hamiltonian 
"Hui(x, t) in Eq. l|23[l are thus replaced by 



q? (x',t f ) = e^ mi(t ^^ x cp V0 (x',t f ) 
Hir>(x' ,t) = Hm (x 1 , t) + mg sin jx' 



(26) 



In parallel with the vertical and oblique average ener- 
gies, Figure also shows the total average energy (E) of 
the trapped atoms after the splitting. This quantity is 
calculated from (E ), (E\), and the total splitting prob- 
ability P s [Eq.(E2)] 



(E) = (1-P S )(E )+P S (E 1 ) 



(27) 



For Ui = 0, no deviation of the cloud is observed, 
and we obtain P s = and (E) = (Eo). This average 
transverse energy is in fact equal to the initial transverse 
energy of the trapped atoms (17.5 /iK). This behavior can 
be seen on the left part of Figure 

When U\ increases by a small amount (U\ < Uq) 
the highest vibrational levels of the vertical guide are 
deviated in the oblique potential (see Figure H3 for in- 
stance), and the average transverse energy of the atoms 



remaining in the vertical guide therefore decreases. A 
striking counter-intuitive effect is that the atoms which 
are deviated also have a translational energy which is 
smaller than the initial average energy of the trapped 
atoms. This happens because these high vibrational lev- 
els are now trapped in a weakly binding potential of 
depth U\ ^ Uq. As a consequence, the total average 
translational energy (E) of the atoms in their transverse 
direction decreases after the splitting of the cloud. With 
the parameters used in Figure El a minimum energy of 
10.9/iK is obtained for XJ\ ~ 0.75 Uo, to be compared 
with the initial average energy of about 17.5 /iK. A sig- 
nificant cooling effect is therefore obtained in the trans- 
verse direction, at the cost of a significant heating in the 
vertical direction. 

Finally, on the right hand side of this figure, with very 
deep oblique potentials (U% ^ 2C/ ), the atom tempera- 
ture in the transverse direction exceeds the initial average 
energy of 17.5 /iK. In this case, a heating process takes 
place due to the fact that the atoms are now trapped in 
a much deeper potential. 



IV. CONCLUSION 

We have proposed a theoretical model for the study of a 
thermal ensemble of cold atoms in a beam splitter device. 
This model has a wide range of possible applications. For 
instance it could be used to describe the dynamics of 
cold atoms trapped and manipulated with the magnetic 
fields created by atom chips. We have used our time- 
dependent semi-classical model to describe the atomic 
dynamics in the presence of two crossing dipole guides. 
We have taken into account the gravity, as well as the 
thermal population of the initial atomic cloud in order 
to compute the splitting efficiency of the beam splitter. 

Our results are in good agreement with experimen- 
tal measurements, and we have presented the influence 
of the main parameters on the atomic dynamics in this 
guiding and splitting configuration. We have shown that 
some eigenstates of the system split coherently in the two 
branches of the guide, and that different average temper- 
atures can be obtained in the different arms of the beam 
splitter. An efficient cooling of the atoms is also predicted 
in the transverse direction. 

All these results indicate that a high degree of control 
can be achieved in this type of cold atom beam split- 
ters, using simple Gaussian laser beams. In the future 
our investigations will concentrate on the theoretical de- 
scription of atom optics devices (guides, mirrors, . . . ) for 
the manipulation of both thermal and coherent sources of 
atoms. Combined with the very impressive capabilities 
of spatial light modulators [U, these techniques should 
effectively allow for the implementation of new exciting 
experimental schemes in the domain of atom optics and 
matter- wave interferometry. 
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